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ABSTRACT 

Author's weighting factor to correct photographic bias due to mass and velocity was 
applied to Smithsonian data for over 2000 sporadic meteors. The factor was supported 
by confirmation of predicted agreement between a velocity histogram for meteors of all 
masses obtained with its use and a raw distribution for large masses. An equation was 
found for geocentric velocity distribution unaccelerated by Earth gravity. Exponent of 
mass of 1. 34 was obtained in equation for cumulative mass distribution using mass 
scale of McCrosky and Posen. An average velocity of 13. 8 km/sec was obtained, vir- 
tually independent of mass. 
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UNACCELERATED GEOCENTRIC VELOCITIES AND INFLUX RATES OF 
SPORADIC PHOTOGRAPHIC METEORS 
byC. D. Miller 
Lewis Research Center 

SUMMARY 


Use was made of an earlier published weighting factor, developed at NASA Lewis 
Research Center, for correction of photographic bias caused by variation in mass and 
velocity of meteors. The weighting factor was applied to all sporadic meteors among 
approximately 2400 for which photographically derived data had been published by the 
Smithsonian Astrophysical Observatory. 

For unaccelerated geocentric velocities of meteors intersecting Earth's orbit the 
results indicated an average value of 13. 80 kilometers per second. This value contrasts 
with previously estimated values ranging from 20 to 40 kilometers per second. Appli- 
cation of the weighting factor showed a close approximation to a log-normal geocentric 
velocity distribution, but with the values of probability density offset 1. 5 kilometers per 
second upward on the velocity scale. As predicted, the weighted velocity distribution 
agreed closely with a raw velocity distribution for a part of the total sample of photo- 
graphic meteors consisting approximately of the heaviest 10 percent. Such prediction 
was made on the assumptions that velocity distribution should be nearly independent of 
mass and that no difficulty should be encountered in photographing even the slowest 
meteors of very large mass. The weighted velocity distribution for all masses was also 
shown to agree with the raw distributions for groups of meteors of lower masses, but 
only for meteors above some critical velocity level at which the slowest of meteors could 
be easily photographed. 

Comparison of semilogarithmic plots of the weighted distribution for all masses with 
the raw distribution for a mass group yielded a correction factor by which the actual 
count of meteors within the mass group could be multiplied to compensate for the failure 
to photograph slower meteors. By use of such correction factors, a corrected distribu- 
tion of mass influx rates was obtained. Such distribution proved to be a straight line on 
a full logari thm ic plot at least down to a mass of 0. 02 gram. The slope of the line and its 
level, perhaps fortuitously, agreed with earlier values within the probable accuracy of the 
method. 

The previously developed weighting factor is believed to be well supported by the 
correlations it produces. 



INTRODUCTION 


A part of the solar system about which existing knowledge is quite limited is a vast 
number of sporadic meteoroids. More extensive knowledge concerning them seems to be 
desirable, not only in keeping with the objective of increasing man’s knowledge of his 
environment, but also because the possibility of collisions with them poses a potentially 
important problem in the planning of long space missions. 

Meteoroids are small particles that are know to be present within the free space of 
the solar system and are generally in elliptic orbits about the sun. When such a particle 
encounters the atmosphere of Earth, it almost always possesses a velocity relative to 
the atmosphere at least as great as the velocity of escape from Earth (about 11 km/sec) 
and may have a velocity above 70 kilometers per second. If the mass of the particle is 
much greater than 10 gram, a trail may be created in the atmosphere sufficiently 
luminous to be visible as a streak across the sky. The appearance of such a luminous 
streak is known as a ’’meteor. ” The particle that causes it, however, is still properly 
termed a ’’meteoroid.” 

-2 

If such a meteoroid has an original mass as great as about 10 gram, the resulting 

meteor may be sufficiently luminous that it can be photographed. The mass range of 
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so-called ’’photographic” meteors extends from this level of about 10" to 10 or 10 
grams. Above the higher levels, meteors may be so rare that a sample of reasonable 
size might contain none, and overexposure of the photographic plate might cause reduc- 
tion to be impossible. 

The mass range covered by the photographic meteoroids may be of particular impor- 
tance for long space missions. A meteoroid having a mass within this range could 
penetrate a substantial thickness of armor. Moreover, the probable frequency of en- 
counter of such meteoroids may be great enough to represent a substantial hazard to a 
space vehicle or space radiator having several hundreds of square feet of vulnerable 
area exposed over a period of many months. 

According to recent results of hypervelocity impact tests (ref. 1), the armor thick- 
ness that can be penetrated by a particle varies as the one-third power of particle mass 
and as the first power of impact velocity. Hence, in assessing the degree of hazard to 
a long-range space mission represented by meteoroids within the photographic range, it 
is necessary to know accurately the influx rate at each possible combination of mass and 
velocity. In other words, it is necessary to know the influx rate for all mass ranges and 
the velocity distribution within each mass range. 

Determination of either influx rates for various mass ranges or velocity distribution 
within a mass range is complicated by a bias due to the effect of either mass or velocity 
on the brilliance of a meteor. Because of this biasing effect, a heavier meteroid is more 
likely to result in a photographic meteor than is a lighter meteoroid, and a meteoroid 
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possessing a higher velocity relative to Earth atmosphere is much more likely to produce 
a photographic meteor than a slower meteoroid. 

A recently published analysis (ref. 2) has shown that for sporadic meteors both the 
mass and velocity biases can be compensated by use of a weighting factor approximately 
proportional to the -4. 22 power of the velocity of a meteoroid relative to Earth atmos- 
phere. This factor is also a function, to a lesser extent, of the angle of meteoroid path 
to the zenith. In use of such weighting factor each photographic meteor is counted, not 
as one meteor, but as one multiplied by the weighting factor as computed from the 
velocity and zenith angle of that meteor. The resulting corrected counts for various 
velocities relative to Earth atmosphere should provide a true representation of velocity 
distribution for all photographic meteors of mass greater than any specified level. The 
weighting factor of reference 2 was first determined theoretically and then revised on an 
empirical basis. The empirical revision, however, agreed closely with the theoretically 
derived factor as to the exponent of atmospheric velocity, which, in practice, is much 
more important than the dependence of the weighting factor on zenith angle. 

An earlier weighting factor that has been widely used for a similar purpose (ref. 3) 
involved atmospheric velocity of a meteoroid only to the negative second power. With 
use of that earlier factor, average velocities of meteoroids relative to Earth atmosphere 
have been estimated at values ranging from 20 (ref. 4) to 40 kilometers per second 
(ref. 5). The much greater negative exponent of atmospheric velocity in the weighting 
factor of reference 2 should be expected to cause a drastic reduction in estimated 
average velocity because it gives much greater weight to the counts of slower particles. 

The magnitude and direction of the effect of the new weighting factor on estimates of 
influx rates of meteoroids of various masses, if any such effect exists, it is not imme- 
diately obvious. However, once the influx rate at various mass levels is determined, a 
calculation of average damage potential would obviously be affected strongly by the 
difference in effective average velocity arising from the difference between -2 and -4. 22 
as an exponent of atmospheric velocity in the weighting factor. 

For this reason, an analysis was undertaken at the NASA Lewis Research Center to 
investigate the effect of the new weighting factor on estimated velocity distribution and 
on estimated influx rates at various mass levels for sporadic meteors. For this 

purpose, the Astrophysical Observatory of Smithsonian Institution provided valuable 
assistance in furnishing on punched cards observational data and various computed 
parameters for approximately 2400 photographic meteors. They were also very helpful 
in providing, in private communication, various items of further information about the 
photographic conditions. The data provided on punched cards were the same as published 
earlier by McCrosky and Posen (ref. 6) . 

Although the weighting factor must be applied to meteor velocity relative to Earth's 
atmosphere, the statistical results of the analysis were sought in terms of geocentric 
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velocities possessed by meteoroids before acceleration by Earth's gravity; that is, an 
overall influx rate and a distribution of influx rate relative to velocity were sought at a 
position in the Earth’s orbit but far from Earth. Such results were desired, instead 
of real conditions near Earth, because (1) a description of the influx distributions should 
be mathematically simpler without the effect of Earth's gravitation and, hence, more 
readily susceptible of expression in the form of simple equations, (2) the influx distri- 
butions without the effect of Earth's gravitation would be encountered for long periods of 
time by a long-range space vehicle, and (3) on the basis of resulting equations repre- 
senting the conditions before acceleration by Earth's gravity the conditions near Earth 
should be analytically derivable without need of further detailed statistical study. 

The analysis was conducted in four stages. In the first stage, a computer program 
used the weighting factor and operated on the data for sporadic meteors of whatever mass 
to yield an adjusted velocity distribution. In the second stage, a similar computer pro- 
gram was used to obtain a velocity distribution for each of 10 intervals of mass, but 
without benefit of the weighting factor. In the third stage of the analysis, the result from 
stage 1 was compared with the result for the heaviest meteoroids from stage 2 to confirm 
or disprove an anticipated agreement. In stage 4 of the analysis, a distribution of me- 
teoroid influx rates relative to mass was derived, with use of a method for correcting 
the count of meteors within each mass group for failures of observation due to insufficient 
brilliance of meteor trails. 

INFLUX DISTRIBUTION RELATIVE TO VELOCITY 

For the first stage of this analysis the theoretical basis will first be described, then 
the computer program will be explained briefly, and finally the resulting velocity 
distribution will be presented and discussed. Although most symbols are defined in 
appendix A, in general, each symbol will also be explained where first used. 


Theory of Adjustment of Velocity Distribution 

The weighting factor that was developed in reference 2 for use in determining adjusted 
velocity distributions is expressed by the equation 

, w(1 ) = (cosz R ,-0.^ F ^a,30 v -4. 22 ft) 
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where Z R is the angle between the meteoroid path and the zenith at the position of first 

encounter of the meteoroid with Earth atmosphere, v ro is the velocity of the meteoroid 

relative to Earth atmosphere at the same position, and F(Z^) av is a statistical average 

function of Z R and meteor position within the field of view of the cameras . The function 

F(Z r ) was designed for the purpose of accounting for the random variation of position 
rt av 

within field of view at which a meteor might appear. It was discussed in detail in refer- 
ence 2 and, hence, will not be further explained here. For all meteors reported in 
reference 6, which formed the basis of the study to be reported here, values of F(Zp) ay 
were computed as part of the work reported in reference 2. The results remained 
available for use in the present study. 

For study of velocity distributions of meteors as they encounter the atmosphere over 
New Mexico, the proper weighting of photographic meteors would be obtained by counting 
each meteor as one impact multiplied by the factor of equation (1). For the work 

reported herein, however, the factor was modified to include an adjustment for gravi- 
tational focussing, that is, the effect of Earth's gravitation in focussing the movement of 
particles in such a way as to increase the influx rate for all particles and to increase the 
influx rates for slower particles in comparison with faster ones. This additional adjust- 
ment was needed because, as explained earlier, it was desired to analyze the conditions 
existing before acceleration of meteoroids by Earth's gravity, rather than the actual 
conditions at Earth's surface. 

The adjustment to compensate for the effect of gravitational focussing was taken from 
an equation developed by Shelton, Stern, and Hale (ref. 7). This equation and the manner 
of its use herein are discussed in appendix B. Its effect is to modify equation (1) to read 


<M2) 


= (cos z R r°- 196 F(z R f- 73 V 4 - 22 



( 2 ) 


where v g is the velocity of escape from Earth at the position where the meteor impacts 
the atmosphere and Vq is the unaccelerated velocity, which the particle would have had 
relative to Earth if Earth's gravity had not affected it. 

Although <P W (2) as determined by equation (2) will provide adjustment of velocity 
distributions to the condition that would have obtained without the influence of Earth 
gravity, it still applies basically to the influx condition as it would then exist in the atmos 
phere over New Mexico, and predominantly during the later night-time hours. A statis- 
tically meaningful correction for the spacewise bias that exists would have required 
rejection of many meteors that were otherwise useable for statistical purposes because 
they came from logical subdivisions of the surrounding space that could not be com- 
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pletely sampled by the atmosphere over New Mexico during night-time hours. Such 
meteors, for example, arrive at the Earth atmosphere from directions south of the 
ecliptic or from directions relatively near that of the sun. It was believed the effect of 
the spacewise bias would not be great. For this reason and because of a belief that any 
attempt to correct the spacewise bias would unduly complicate this pr elim i na ry work, 
a decision was made to ignore the spacewise bias involved in the restriction of the 
exposure area to New Mexico during the later night-time hours. 

A third weighting factor according to the equation 


^w(3) - 



( 3 ) 


was used for obtaining velocity distributions adjusted for gravitational focussing but not 
for the biasing effect of mass or the velocity v . For convenience in the following 
description of the general use of a computer program, a unity weighting factor is defined 
as 


^w(4) = 1 W 

A computer program was written for the preparation of tabulations of values of Vq 
with use of the weighting factor <P w ( n )> with a value of 2, 3, or 4 for n. Data cards fur- 
nished by Smithsonian Institution Astrophvsical Observatory were read and re-recorded 
for use with the computer program in a manner to ensure that the data for any pertinent 
meteor could not be overlooked on any run of the program. The cards contained exactly 
the data published in reference 6. They were checked with reference 6 for completeness. 
Data cards were excluded for meteors for which not all the needed data had been deter- 
mined. In execution of the computer program, the following operations were performed 
for a single meteor, then for the next meteor, and so on. 

(1) The parameters , Vq, and cos Z R’ were read, as well as serial number and 
the notation of shower association if present. (The parameters v m , Vq, and cos Z R 
correspond, respectively, to the column headings v w , Vq, and CZ R in table 7 of 
reference 6.) The value of F(Z R ) ay as determined in the work reported in reference 2 
was also read. 

(2) The meteor was rejected from consideration if a shower association was indicated. 

(3) The meteor was rejected from consideration if the value of cos Z R was less than 
0.20. (It was not believed this limitation would be statistically important; it resulted in 
rejection of approximately 20 meteors, or less than 1 percent of the total. This provision 
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was equivalent to limiting the acceptance of the atmosphere to a cone of approximately 
102. 5° included angle rather than 180°. The restriction was necessary because of the 
sensitivity of the weighting factor cp or <p w ( 2 ) to small errors in the smaller values 
of cos Zj^. In particular, these factors go to infinity for cos Z R = 0. 00. As several 
meteors were on the list with cos Z R = 0. 00, they would have made use of these weight- 
ing factors quite impossible. ) 

(4) The value of was computed with a value of 11. 1 kilometers per second for 

v_ (at an altitude of 90 km). 

(5) The integer next higher than Vq in kilometers per second was found and the 

value of was added to the content of a storage location corresponding to that 

integer (initialized to zero at the start of the program). 

(6) The value of was multiplied by the value of Vq and added to the content 

of a storage location (initialized at zero) for later use in determination of average 
velocity. 

After these operations had been performed for all meteors, the content of the storage 
location for each integral value of Vq was divided by the sum of the contents of all and 
multiplied by 100 to give the adjusted percentage of total count for each velocity interval. 
(Hereafter, this procedure for reduction of total counts to unity, or to 100 percent, will 
be referred to as "normalization. " That term will also be used relative to the provi- 
sion for unity, or 100 percent, area under a distribution curve. Factors or divisors 
used for the purpose will be referred to as normalizing factors or divisors. ) The content 
of the storage location for use in determination of average velocity was also divided by 
the sum of the contents of all storage locations for integral values of v^, to yield the 
average value of velocity. 


General Nature of Velocity Distribution 

Results from runs of the program that has just been described, with n equal to 2 and 
with n equal to 4, are presented in figures 1 and 2 with several different plotting scales. 
These plots are actually histograms. However, for these figures and others, straight 
lines are drawn from one value to the next rather than the customary horizontal lines 
forming the tops of the bars. This change was made to avoid confusion in the comparison 
of two histograms in the same figure. 

The raw distribution in the linear plot in figure 1(a) represents results obtained 
from the computer program with use of * or eac h meteor, that is, with no ad- 

justment even for gravitational focussing. It shows the familiar bimodal distribution with 
two peaks: one for direct orbits in the neighborhood of 14 kilometers per second, the 
other for retrograde orbits in the neighborhood of 60 to 70 kilometers per second. This 
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0 10 20 30 40 50 60 70 80 

Unaccelerated geocentric velocity, vq, km/sec 

(a) Linear scale. 



Unaccelerated geocentric velocity, vq, km/sec 


lb) Semi logarithmic scale. 


Figure 1. - Comparison of raw distribution of velocity of meteors with distribution adjusted 
for photographic bias of meteoroid mass and velocity and for gravitational focussing. 
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Geocentric velocity, v G , km/sec 
(Offset +1.5 km/sec) 

1 1 i i i i j 

1 2 4 10 20 40 100 

True logarithmic scale 

Figure 2. - Fit of log-normal equation to velocity distribu- 
tion of meteors corrected for gravitational focussing and 
bias due to meteoroid mass and velocity. 


distribution superficially indicates that particles in retrograde orbits might well repre- 
sent the principal hazard because of their high velocities and substantial frequency. 

However, the adjusted distribution in the same part of figure 1, obtained with use of 
<p w ( 2 )> shows that the retrograde orbits are comparatively rare and that the frequencies 
of velocities in the neighborhood of 8 to 18 kilometers per second are very high in com- 
parison with all other velocities. The run of the program showed an average adjusted 
velocity (with <p w ( 2 )) °1 13. 8 kilometers per second as compared with a raw average 
(with of 30. 0 kilometers per second. 

The same results as plotted linearly in figure 1(a) are plotted on a similogarithmic 
scale in figure 1(b). Here, it is seen that the peak of the adjusted distribution for 
retrograde orbits is lower by more than two orders of magnitude than the peak for the 
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direct orbits. Although the adjusted distribution is still bimodal, the departure from a 
single mode is not nearly so great as in the raw distribution. It is seen that both raw 
and adjusted curves drop to very low levels in the neighborhood of 1 to 2 kilometers per 
second. Such particles would have to have orbits almost the same as that of the Earth 
and nearly the same velocity relative to the sun. It will be seen that, below 3 kilometers 
per second, the adjusted distribution falls under the raw distribution. This condition 
exists because, in this velocity range, the adjustment for gravitational focussing that 
is included within the function <P W ( 2 ) exceeds the adjustment for effect of mass and 
velocity on luminosity. 

The same results for adjusted velocity distribution as shown in figure 1 are plotted 
in figure 2 on a full logarithmic scale. The values shown on the abscissa scale are true 
values of the midpoints of the 1-kilometer-per-second velocity intervals shown on the 
abscissa scales in figure 1. However, they are shifted 1. 5 kilometers per second toward 
the right relative to the logarithmic scale. This shift resulted from a trial -and-error 
procedure to obtain the best possible fit of a curve representing a distribution equation to 
the plotted points. The results, as delivered by the computer and as plotted in figure 1, 
were referred to the maximum of each velocity interval. Hence, the offset of points by 
1. 5 kilometer per second (fig. 2) is an offset of only 1 kilometer per second relative to 
the abscissa scales of figure 1. In cases where two or more data points were too close 
together to plot separately in figure 2, a single point was plotted and the numeral 2 or 3 
was entered in the figure to indicate that such a number of data points correspond to the 
plotted position. 

The curve in figure 2, which appears to represent the plotted results closely, is a 
true parabola. As shown in appendix C, it represents the following offset log-normal 
equation. 


f(v Q ) = 0. 07 exp 



G 


+ 1.5) - log e 11.7 
0. 44 



( 5 ) 


So it is seen that equation (5), plotted as the parabola in figure 2, represents a tentative 
analytical velocity distribution that closely describes the data. No theoretical reason is 
offered as to why the distribution should be log-normal with an offset of 1. 5 kilometer per 
second; the discovery is empirical. 
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For geocentric velocities below 30 kilometers per second, equation (5) represents 
very roughly the same normal distribution as reported by Wall (ref. 9), which he derived 
by other methods. At higher velocities, the distribution according to equation (5) differs 
greatly from his. 1 


RAW VELOCITY DISTRIBUTIONS FOR MASS GROUPS 

In the second stage of the analysis, unaccelerated velocity distributions (histograms) 
were obtained for each of 10 mass groups with use of weighting factor cp^^y 111686 
histograms will be used for comparisons with the histogram for all sporadic meteors, 
adjusted with weighting factor (? w ( 2 )> b y which the most essential objectives of this 
analysis will be accomplished. Hereinafter, the distributions for the 10 mass groups, 
with adjustment only for gravitational focussing, will be referred to for brevity as 
’’velocity distribution for mass group 1, ” ’’velocity distribution for mass group 2, ” and 
so on, or collectively ’’velocity distributions for mass groups. ” The distribution earlier 
described, obtained for sporadic meteors of all masses with use of weighting factor 
<p w ( 2 ) (plotted as the square symbols in fig. 1 and as the round symbols in fig. 2), will 
be referred to hereinafter as the ’’velocity distribution for all masses. ” 

To obtain the velocity distributions for mass groups, the data cards for all photo- 
graphed meteors were sorted to eliminate cards for which mass was not given and to 
arrange the remaining cards in order of descending mass. The cards were then phys- 
ically divided into 10 mass groups, each containing the same number of cards ±1/2 card. 
The computer program earlier described was used to tabulate the distribution of velocity 
Vq for each of the 10 mass groups with use of weighting factor (P-^^y excluding each 
shower -associated meteor as before. The percentage of shower -associated meteors 
proved to be substantially higher for low masses than high, so that the meteors left in 
each mass group after rejection of the shower -associated meteors ranged from 217 in 
the heaviest group to 172 in the lightest. However, the lighter groups were still thought 
to be large enough to be statistically significant. The groups were used without change, 
in preference to redividing them with a consequent extension of the higher mass ranges 
downward into the lower masses. 

The velocity distribution for each of the 10 mass groups resulting from the computer 
tabulations is plotted in figure 3. Also plotted in figure 3 is the velocity distribution for 
all masses. 

The velocity distributions for the 10 mass groups are normalized with implicit use 
of a divisor believed to be equal to the true total influx rates for those groups, which in 
most cases differed from the observed influx rates. The manner of this normalization 


11 


Is 




Percentage of total count having velocity equal to abscissa, _ j km/sec 
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Figure 3. - Continued. 
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(g) Mass group 7. 
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(h) Mass group 8. 

Figure 3. - Continued. 
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Unaccelerated geocentric velocity, v G , km/sec 
(1) Mass group 9. 
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Unaccelerated geocentric velocity, v G , km/sec 


<j) Mass group 10. 
Figure 3. - Concluded. 



will be described somewhat later. In the single case of mass group 5 (fig. 3(e)), the 
distribution for the group is also shown normalized by using a divisor equal to the ob- 
served influx rate corrected only for the effect of gravitational focussing. The additional 
plotting is shown in this part of figure 3 because it was chosen as the most convenient 
figure for later detailed discussion. 

The velocity distributions for the mass groups in general cannot be regarded as 
realistic, because they have no correction for the velocity bias that existed in the photo- 
graphic process. Their only purpose is for use in the following analytical comparisons. 


RAW VELOCITY DISTRIBUTION FOR HEAVIEST METEOROIDS 

The third stage of the analysis involves investigation of the prediction that an ad- 
justed velocity distribution should not differ from an unadjusted distribution for a very 
heavy mass group whose lower limit of mass will be termed m m j n - In reference 2, in 
the development of the expression for > an assumption was used to the effect that 

the velocity distribution is substantially the same for any mass within the photographic 
range. Such assumption appears to be reasonable (1) because the mass of a particle 
does not enter into the equations that interrelate orbital parameters (including velocity) 
and (2) because no theoretical mechanism is known to produce substantial orbital changes 
on a mass basis within the mass range involved in this study, except possibly at very low 
geocentric velocities (circular orbits) where the Poynting- Robertson effect (ref. 10) might 
assume some importance. 

If the assumption of the same velocity distribution for any mass is correct and if 
equation (1) is also correct, it should be expected that the velocity distribution for all 
masses, determined with use of weighting factor m ight be identical with a velocity 

distribution for a very heavy mass group, obtained with use of weighting factor (p w ^y 
For, if the value m m j n within the group were great enough, virtually all meteoroids 
within the group should produce trails that could be easily photographed regardless of 
velocity. The minimum practical limit of velocity of impact on the atmosphere should be 
the velocity of escape from Earth. Hence, if the value of m m £ n were great enough to 
allow easy photography at the velocity of escape, few meteoroids at any impact velocity 
would escape detection. 

Striking confirmation of the prediction described is shown in a comparison of results 
presented in figure 4. The dashed-line histogram in that figure represents the average 
distribution of velocities Vq for the two heaviest mass groups, normalized relative to 
the observed count of meteors adjusted with use of weighting factor <p w Qy Initially, the 
decision was made to use the total content of the first two mass groups for the compar- 
ison in this figure because the histograms for these two groups looked very much alike. 
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Figure 4. - Comparison of distribution of velocity of meteors within combination of two 
heaviest mass groups with distribution for all masses adjusted for photographic bias due 
to meteoroid mass and velocity. (Both plots are adjusted for effect of gravitational focussing). 


A more elaborate indication that this choice was a good one will be mentioned later. The 
solid-line histogram in figure 4 is the velocity distribution for all masses, as earlier 
defined and as sketched earlier in figure 1(b). 

Throughout the velocity range from about 20 to 79 kilometers per second, the two 
histograms show the same general trend, but with the histogram for the two heaviest 
mass groups generally on a slightly higher level. In the velocity range from 1 to about 
10 kilometers per second the trends are again similar, but with the relative levels re- 
versed. The differences of level are slight, and the comparison, therefore, amounts to 
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a general confirmation of the validity of the adjusting factor • Specifically, the 

appearance of the two histograms in figure 4 for geocentric velocities above 20 kilometers 
per second is a confirmation, as exact as the point-to-point irregularities permit, of 
the assumption used in reference 2 that the same velocity distribution exists for all 
masses. 

Two possible causes of the failure of the adjusting factor <P W (2) to P rovide exact 
agreement between the two histograms will be discussed. 

(1) Some failure to observe low -velocity meteors because of faintness of the traces 
in the photographic plates may have occurred even within the heaviest mass groups. If 
the meteor sample had been larger, the heaviest groups could have been further subdivided 
to determine whether the raw histograms for the heaviest subdivision might have dupli- 
cated the adjusted histogram for all masses. But such subdivisions would not have pro- 
vided a large enough sample within each subdivision to serve as an adequate statistical 
base. One slight indication contrary to this explanation is the close agreement between 
the raw histograms for mass groups 1 and 2 in figures 3(a) and (b). 

(2) In the derivation of the factor in reference 2, it was not possible to take 

into account the effect of meteoroid fragmentation. As discussed in that reference, 
fragmentation should tend to increase the brilliance of a meteor. Fragmentation is also 
known generally to cause a more pronounced effect with meteors of lower speed (refs. 11 
and 12). The factors and ¥^(2) were designed to magnify real counts of meteors 

in a manner dependent, among other things, on their faintness. If, due to fragmentation, 
the meteors of lower speeds are not actually as faint as they were treated in the derivation 
of the factor > then that factor would excessively magnify the count of meteors in 

the low-speed range and cause the slight systematic discrepancy between the two histo- 
grams that is seen in figure 4. 

The question of disagreement of these histograms in the low-velocity region can not 
be resolved here, but must await later study. The assumption of a single velocity dis- 
tribution for all masses and the correctness of the weighting factor <p w ^ for velocities 
lower than about 20 kilometers per second are only approximately confirmed. But they 
are not contradicted by the appearance of the histograms, even on an exact basis. 

Hereinafter explanation (2) for this anomaly will be treated as the correct expla- 
nation, but with the reservation that later research may show explanation (1) to be cor- 
rect. No significant change in final results would then be required. 

INFLUX DISTRIBUTION RELATIVE TO MASS 

Determination of an influx distribution relative to mass will now be undertaken as 
the last stage of this analysis. Hereinafter, that influx distribution will be called "mass 
distribution. " 
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Meteoroid Mass Determination 


In determination of a mass distribution from the data of reference 6, some uncer- 
tainty exists concerning a basic reference level. The reference level used here places 
the upper and lower limits of the 10 mass groups at values that do not relate directly 
to the masses given in reference 6. They have been adjusted upward from the values of 
reference 6 by a factor 6.4. 

Since the publication of reference 6, the most probable value of luminosity coefficient 
of meteors has been revised. As stated in that reference, meteoroid mass m was 
determined by integration of the equation 


m = 



( 6 ) 


where Tq is the luminosity coefficient used in reference 6, v is the instantaneous 
velocity, T is the lifetime of the meteor, and I is the instantaneous luminous intensity. 
It is apparent, therefore, that a corrected mass m cor should be expressed as 


m cor = m 


T 0(cor) 


( 7 ) 


where To(cor) * s the value of luminosity coefficient considered to be most probable at 
this time. 

From information given in references 6 and 13 to 15, a value 


— T ° — = 6.4 (8) 

r 0(cor) 

was deduced for use in equation (7) for the purposes of the present analysis. After com- 
pletion of the analysis, but before publication, reference 16 became available. On page 4 
of that reference there is indication that the ratio should be 

T ° =6.46 (9) 

r 0(cor) 


20 



As equations (8) and (9) are nearly identical, further justification of the factor 6.4 used 
in the present analysis will not be included herein. 


Correction of Counts of Meteors 

A method was found for correction of the counts of meteors within the 10 mass 
groups for failure of observation. The method involves a comparison between the 
velocity distribution for all masses as earlier defined (the solid- line histogram of figs. 1(b) 
and 4) and each of the 10 dashed-line histograms (fig. 3) . Those 10 histograms, as 
explained earlier, represent the unaccelerated velocity distribution for each of the 
10 mass groups, adjusted only for gravitational focussing. 

In the comparison a choice was necessary, whether to use the solid-line histogram 
of figures 1(b) and 4 or the dashed-line histogram of figure 4. The decision was made to 
use the histogram first mentioned because (1) only the parts of the histogram above 
about 20 kilometers per second will have importance in the comparisons and (2) in all 
velocity regions above about 20 kilometers per second, the histogram used for compari- 
son needs to be as accurate as possible. The dashed-line histogram in figure 4 is not 
dependable above about 35 kilometers per second. Above the geocentric velocity of 35 
kilometers per second, in the two groups of greatest mass, meteors were found only at 
velocities of 36, 38, 45, 48, 55, and 69 kilometers per second. A solution of the prob- 
lem of plotting the percentages of total influx at these isolated velocities in such a 
manner that the values would be significant in relation to the other parts of the histogram 
involves unavoidable compromise and uncertainty. (This problem also exists at both 
high and low velocities for the dashed-line histograms of the various parts of figure 3. ) 

The problem was resolved in an arbitrary and not entirely satisfying manner by assum- 
ing, after normalization, that an isolated count should be distributed impartially between 
the velocity class within which it actually occurred and the velocity classes separating 
it from the nearest nonvacant class in the direction toward the main body of the data. A 
point was actually plotted, however, only for the velocity class within which the count 
occurred. 

The count of meteors within a mass group, before correction, was the total count 
for all velocity classes delivered by the computer program with use of the weighting 
factor <p w (3). This factor provided an adjustment only for gravitational focussing, by 
counting each sporadic meteor as one multiplied by the factor <p w @y The first step in 
the use of the counts for the mass groups was the determination of a basic factor a 1 
by which the counts for the various masses could be multiplied to convert them to par- 
ticles per square meter per second. For this purpose all meteors of any mass of mag- 
nitude equal to or less than zero were counted, without adjustment even for effect of 
gravitational focussing. The resulting total count (436) was used in conjunction with 
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equation (16) of reference 17 which reads as follows: . 

lo SlO N 0 = _4 * 33 (1°) 

where N Q is frequency of occurrence of meteors of magnitude zero or brighter per 
square kilometer per hour. Such equation was accepted as correct for the purpose of 
this analysis. On such basis, a value of or* by which any count may be multiplied to 
convert it to particles per square meter per second is 

a* = & 6S<10 ^ * 10 6 = 2. 982XKT 7 (11) 

436 x 3600 


that is, 

C = C V (12) 

where C T is a count of meteors of a given mass range adjusted for gravity focussing or 
not and C is a corresponding frequency of influx in particles per square meter per 
second. 

As mentioned earlier, no adjustment for spacewise bias in the photography was made 
for the present analysis. It appears, also, that no adjustment for the spacewise bias: 
was made in the work of reference 17. Hence, equation (11) will be applicable to the 
atmosphere over New Mexico during the late night-time hours rather than to a randomly 
oriented surface in space. The equation will therefore be in harmony with the velocity 
distributions and the counts of meteors of various masses as they will be used in the 
present analysis. 

In the description of the next steps in the procedure, involving adjustment of the 
counts, the fifth mass group will be used as an illustrative example. It was chosen for 
discussion because the features to be studied are more clear cut than in other cases. 
However, exactly the same procedure as that which will be described relative to mass 
group 5 was used for each of the mass groups in turn. For each mass group i, the 
computer count of meteors adjusted with weighting factor <P W ( 3 ) will be designated 

C a(i)' 

The solid-line histogram in figure 3(e) shows the velocity distribution for all 
masses. It is the same as the solid-line histogram of figure 1(b) or figure 4. 

The shape of the velocity distribution for mass group 5, with adjustment for gravity 
focussing, appears as two dashed-line histograms in figure 3(e), on levels corresponding 
to the use of two different normalizing divisors as earlier described. These two 
differently normalized histograms will be discussed in more detail. 
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Before the construction of figure 3(e) the normalizing divisor for the triangular 
points (132. 4) was known, simply as the computer’s count of meteors adjusted only for 
gravitational focussing. The normalizing factor for the circular points was determined 
by the construction of the figure which involved the implicit use, mentioned earlier, of 
a normalizing factor believed equal to the true influx for the mass group. 

If the histogram for all masses (the square points) is actually correct as the veloc- 
ity distribution for the fifth mass group, then it should agree with the histogram for 
mass group 5 as to shape on the semilogarithmic plot everywhere to the right of some 
critical geocentric velocity. That critical geocentric velocity should be the velocity at 
which even the particle of smallest mass within the group would be bright enough to be 
discovered without failure on the photographic plates. At lower velocities, at least 
occasional failure would occur in the discovery of meteors even though they might have 
great enough mass to place them within the group. At those lower velocities, the failure 
to discover faint meteors would cause the histogram for mass group 5 to be on too low a 
level in comparison with the higher velocities. The histogram would therefore not have 
the same shape as the solid-line histogram at these lower velocities. 

In consequence, the correct level for the histogram representing the velocity dis- 
tribution for mass group 5 may be found simply by moving the histogram downward 
without change of shape from the original position represented by the triangular symbols 
to the position shown by the circular symbols. The correct level should be reached when 
the maximum agreement exists between the solid-line histogram and the dashed-line 
histogram within the parts of those histograms generally toward the right. Such a 
maximum agreement was thought to exist when the dashed-line histogram reached the 
level shown by the circular symbols, on a lower level than that shown by the triangular 
symbols by a distance equal to log 6.70. 

At the lower level, the velocity distribution for all masses (square points) coincides 
with that for mass group 5 (circular points), in general level and in shape, as exactly as 
could be desired throughout the velocity range from about 28 to about 58 kilometers per 
second. Below that velocity range, the histogram for mass group 5 falls off below that 
for all masses because of the marginal luminosity for particles in this mass group at 
these lower velocities. On the other hand, at velocities above about 39 kilometers per 
second the population density of this mass group was so small (10 meteors total) that 
great weight could not have been given to the level of the dashed-line histogram within 
that velocity range, even if it had not agreed well with the solid-line histogram. 

From inspection of figure 3(e) the critical velocity above which the histograms agree, 
for the fifth group, is about 28 kilometers per second. As demonstrated in appendix D, 
the distance of the downward shift of the dashed-line histogram necessary to bring about 
agreement of the two histograms throughout the velocity region above the critical should 
be the logarithm of a correcting factor cp by which the count of particles in the mass 
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group should be multiplied to make that count correct; that is, with substitution of 
for <p of appendix D, for mass group i, 

T f o i 

C? = C a ^a <p^ particles m sec - (13) 

where C? is an influx rate for mass group i corrected for the shift of the histogram . 
The amount of the necessary shift in the case of mass group 5 was measured as 

log ^( 5 ) = log 6. 70 (14) 

while the value of C a ^ delivered by the computer was 132. 4 meteors. Consequently, 
by equation (13) the count for mass group 5 corrected for the shift of the histogram 
should be 

Cg = 132.4 <*V( 5 ) = 2.65x10”*^ particles m - ^ sec - * (15) 

All other parts of figure 3 show the same solid-line histogram as in figures 1(b), 4, 
and 3(e). Also, each part of figure 3 shows a dashed-line histogram for one of the mass 
groups from 1 to 10. In each case, as has been explained in reference to figure 3(e), the 
dashed-line histogram has been shifted downward by the necessary amount to bring it 
into agreement with the solid-line histogram throughout the velocity range above which 
particles are not lost because of faintness of exposure. In the downward shifting of the 
dashed-line histograms, little weight was given to the highest velocity regions if the 
population density in such regions was sparce. In each case, as with figure 3(e), the 
distance of the downward shift is the logarithm of a correcting factor cp'^ for application 
to the observed count C & ^, for use in equation (13). 

Values of log <p'^ obtained by the method described for the 10 mass groups are 
tabulated in table I. Also shown is the lower limit of mass for each group. 

The count Cg of equation (15) and analogous counts c| for other mass groups, 
corrected for the shifts of the histograms, need an additional correction. From figure 4 
an inference may be drawn that the velocity distribution for all masses is on too low a 
level throughout the velocity range above about 20 kilometers per second, because it 
appears to be on too high a level at lower velocities. Figure 4 should not be interpreted 
as an indication that meteoroids of geocentric velocity less than 10 kilometers per 
second are more rare for greater masses than for lower masses. The meteors of 
these lower velocities that were actually observed were all of large mass. A greater 
actual count at these lower velocities would have raised the level of both histograms in 
figure 4 proportionately, and a lower actual count would have lowered both propor- 
tionately. So instead, as discussed earlier, the condition under discussion will be treated 
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TABLE I. - HISTOGRAM SHIFTS FOR 
VARIOUS MASS GROUPS 



Mass 

group 

Shift required to cause 
histograms to superimpose 
in regions of higher velocity, 
Log 

Lower 

mass 

limit, 

g 

Figure 

1 

log (1. 54) 

0.705 

3(a) 

2 

log (1. 25) 

.397 

3(b) 

3 

log (2.90) 

.250 

3(c) 

4 

log (4. 45) 

. 1603 

3(d) 

5 

log (6.70) 

. 1025 

3(e) 

6 

log (17. 0) 

.0641 

3(f) 

7 

log (29. 0) 

. 03786 

3(g) 

8 

log (72. 0) 

. 01922 

3(h) 

9 

log (150. 0) 

.0077 

3(i) 

1° | 

log (1700) 

. 000635 

3(j) 


as an indication that because of the effect of fragmentation factor or <^(2) no * 

quite correct, at least for comparison of the velocity region below about 10 kilometers 
per second with that above about 20 kilometers per second. The apparent inaccuracy 
causes the adjusted velocity distribution for all masses to show too high an influx rate at 
velocities less than about 10 kilometers per second. Conversely, then, the indicated 
influx rate at velocities higher than about 20 kilometers per second must be too low. The 
effect of this incorrect level on the counts of meteors should now be considered. 

In the preceding parts of the analysis, the incorrect shape of the histogram for all 
masses within the velocity range below 20 kilometers per second has been of no conse- 
quence. For, on examination of appendix D, it may be seen that the entire correct 
histogram for all masses is not needed. The derivation is equally valid if it is assumed 
that only the portion being matched has the correct shape and level. However, because 
of the incorrect level of the utilized part of the histogram for all masses, all the shifts 
of dashed-line histograms that were made in construction of figure 3 were actually 
greater than they should have been by a constant amount. That amount may be estimated 
from the values of log tabulated in table I as follows. 

It is observed that the values of log cp ^ and log (p^y namely , log (1. 54) and log 
(1.25), progress in the opposite sense to the progression from log cp ^ through log 
<p^l 0 ) . The comparison of log and log < pfy would seem to indicate a need for a 
greater adjustment for meteors too faint to photograph in mass group 1 than in mass 
group 2. This condition cannot logically exist. For that reason, mass groups 1 and 2, 
and only they, are believed to include only masses so great that the meteoroids practi- 
cally never escape detection however low their geocentric velocities. (This condition is 
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the earlier mentioned independent confirmation that the total content of mass groups 1 
and 2 was a good choice for use in the comparison that was made in fig. 4.) Accord- 
ingly, log (p^ and log should be equal. The difference between them must be 
due to random irregularities or inexactness of the method. Moreover, if the histogram 
for all masses were correct at geocentric velocities above about 20 kilometers per 
second, both (p'^ and <pfy should equal zero. Hence, the average of these two values 
should be at least a fair measure of the excess shift of the dashed-line histograms, not 
only for groups 1 and 2 but for all 10 mass groups. Accordingly, a corrected shift, 
<P(iy is defined by the equation 


log <p ^ = log <p j.) 


log(l. 54) + log(l. 25) 
2 


(16) 


or 


log <?(.) = log - log(l. 39) (17) 

For any mass group i, the shift together with equation (17) permits the writing 
of an expression for a fully corrected count of meteors for the group as 

C i = C a(i) a,< P(i) (18) 

Thus, in the final result, equation (18) replaces equation (13), C. replaces C? of 
equations (13) and (15), and replaces <pfo in equations (13) to (15). In particular, 
for mass group 5, equation (15) is changed to 

C 5 = C a(5) <*V( 5 ) = 132. 4x2. 982X10" 17 X JliZ® 


= 1. 902x10“^ particles m"^ sec' 


(19) 


As mass group 1 contains all meteors of mass greater than its lower limit, a fully 
corrected cumulative count for any mass group i may be determined with use of 
equation (18) as 


'cum(i) 


• ■£ * 

j=i 


( 20 ) 
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or 


^cum(i) ^cum(i-l) + 

As the same procedure that has been described with mass group 5 as an example 
was used with each mass group 1 to 10 in turn, a value of C cum (i) was determined with 
use of equation (20) for i = 1 and with use of equation (21) thereafter. For mass group 
4 the result was 


Cp um (4) = 2.916X10 -14 particles m -2 sec - * (22) 

Accordingly, for the fifth mass group, from equations (19), (21), and (22) 

C cum(5) = C cum(4) + C 5 = 4 * 81 8>< 10 " 14 particles m -2 sec -1 (23) 


RESULTS 


Tabulated in table II, for each of the 10 mass groups, are value of count of meteors 
C a( i) corrected only for the effect of gravitational focussing, logarithmic correction 


TABLE H. - UNACCELERATED METEOROID INFLUX DISTRIBUTION RELATIVE TO MASS 


Mass 

group 

Count 
corrected 
for gravity 
focussing, 

C a(i) 

Log of 
correction 
factor 

for unobserved 
faint meteors, 
log <p (i) 

Fully 

corrected 

count, 

c i> 

particles 
m” 2 sec" 1 

Fully corrected 
cumulative 
count, 

^cum^)’ 

particles 

m“ 2 sec" 1 

Lower 

mass 

limit, 

g 

Figure 

1 

132.1 

log (1. 11) 

3.94X10' 15 

3.94X10” 15 

0.705 

3(a) 

2 

109.8 

log (0. 90) 

3.28X10" 15 

7.22X10' 15 

.397 

3(b) 

3 

129.6 

log (2. 08) 

8. 06xl0” 15 

1. 528X10” 14 

.250 

3(c) 

4 

145.5 

log (3.20) 

1. 388xl0” 14 

2.916X10” 14 

. 1603 

3(d) 

5 

132.4 

log (4. 82) 

1.902X10" 14 

4. 818X10” 14 

. 1025 

3(e) 

6 

133.2 

log (12.22) 

4. 86xl0” 14 

9.68X10” 14 

.0641 

3(f) 

7 

152.6 

log (20. 85) 

9.48X10" 14 

1.916x10" 13 

. 03786 

3fe) 

8 

139.2 

log (51.8) 

2. 15x10” 13 

4. 07x10” 13 

. 01922 

3(h) 

9 

132.8 

log (108. 0) 

4. 28X10” 13 

8. 35x10” 13 

.0077 

3(i) 

10 

119.7 

log (1223) 

4. 36x10” 12 

5. 20x10” 12 

. 000635 

30) 
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factor <p^ according to equation (17), fully corrected count CL according to equa- 
tion (18), fully corrected cumulative count C cum (j\ according to equation (20) or (21), 
and lower mass limit. In the computation of fully corrected counts C ^ as shown in 
table I, two minor departures were made from the procedure that has been described. 
For mass group 1, the excess value of <p^ above unity is contrary to the earlier con- 
clusion that no particles in mass group 1 could have escaped detection because of low 
velocity. The excess was considered to be due to inaccuracy of the method and the 
value was treated as exactly unity. For mass group 2, the deficit of <p^ below unity 
would correspond to a negative number of particles escaping detection because of low 
velocity. Hence, the deficit was considered to be due to inaccuracy of the method and 
the value was treated as exactly unity . 

The values of fully corrected cumulative count C cum ^ as shown in table II are 
plotted on a full logarithmic scale relative to the group lower limits of particle mass as 
the circular points in figure 5. The straight line in the figure representing the results 
of the present analysis was constructed on the basis of the points for mass groups 1 to 
8. It virtually passes through each of these points. The plotted points for mass groups 
9 and 10 were given no weight in the construction of the line because it was fully expected 
that these points would fall below their correct levels. From examination of figure 3 it 
is clear that the photographic exposure threshold is not sharply defined; that is, instead 
of a sharp cut off below a particular velocity for a given mass, there is a gradual failure 
to detect particles of progressively lower velocity. It follows reasonably that, if all 
particles of mass equal to or greater than the minimum of the 10th mass group had been 



Figure 5. - Distribution of influx rate of meteors relative 
to mass. 
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detected, even at the highest velocities of 77 to 79 kilometers per second, then some 
particles of even smaller mass would have been detected at those velocities. As none 
such were detected, it follows that particles with mass that would have placed them in 
the 10th group and velocities in the neighborhood of 77 to 79 kilometers per second 
must have gone undetected. In comparing the two histograms of figure 3(j), nothing 
better could be done than to place the points of the two histograms level at velocities of 
77, 78, and 79 kilometers per second. Undoubtedly, the points for mass group 10 should 
have gone lower even in that velocity range. Hence, in view of the demonstrated 
absence of a sharp cutoff at the exposure threshold, the 10th point as plotted in figure 5 
should fall below its correct level, and it should not be surprising that the ninth point 
would also do so. 

The slope and position of the line drawn through the circular symbols for the first 
8 mass groups call for a = 10 -1 ^' and 0 = 1. 34, as nearly as can be estimated, in 
the mass distribution equation 


F > = am~P (24) 

In that equation is influx rate of meteors of mass greater than m in particles per 
square meter per second. The value of 0 equal to 1.34 has previously beed widely, 
though not universally, accepted. It is the value for the "Whipple 1963A” line (ref. 13) 
which is also shown in figure 5 for comparison. The close agreement, however, seems 
to be fortuitous. This same value was obtained by Hawkins and Upton (ref. 17), from 
whose work Whipple adopted the slope used for his 1963A line. But Hawkins and Upton 
used photographic magnitudes determined by Hawkins and Southworth (ref. 14) . Kresak 
(ref. 18) has shown an appreciable systematic difference in the magnitude and mass 
scales obtained from references 6 and 14. As the mass values used in the present 
analysis were those of McCrosky and Posen (ref. 6), different results for the parameter 
/3 should be expected between the present analysis and that of Hawkins and Upton. With 
reference to existing dissent relative to the value 1. 34, a value of /3 as low as 0. 88 
has been advanced by Dohnanyi (ref. 19). 

The present analysis makes clear the fact that an accurate determination of the 
mass distribution throughout the photographic range is possible only with explicit use 
of, or implicit determination of an accurate weighting factor to correct the photographic 
bias due to differences of mass and velocity. Prior to reference 2, however, no exten- 
sive analysis appears to have been published in relation to any weighting factor to be 
used for this purpose. Although at least to some extent fortuitously, the value of /3 
equal to 1.34 has been well confirmed by the present analysis, with use of the weighting 
factor from reference 2. Also, the straightness of the cumulative mass distribution 
curve on a full logarithmic plot has largely been demonstrated, whereas such straight- 
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ness had previously been effectively assumed except for the larger masses within the 
photographic range. 

The use of the mass scale of McCrosky and Posen (ref. 6) herein and in reference 2 
does not necessarily involve a preference for this scale relative to that of reference 14. 
This choice was only a matter of convenience. The weighting factor derived in refer- 
ence 2 would be supported by the present analysis, even if the mass used here were not 
correct, so long as the mass used was a monotonic increasing function of the true mass 
alone. If the mass used here were a power of the true mass times a constant, the value 
of /3 derived here could be readily adjusted. 

14 1ft 

The value of or for the ’’Whipple 1963A” line is 10 ' , somewhat higher than 

the value 10 -1 ^'^ from this analysis. The lower level from this analysis is believed to 
be roughly the same as might be expected because of the correction for the effect of 
gravitational focussing. The Whipple line is without an adjustment for gravitational 
focussing. Also, for the line shown, the value of a used by Whipple has been multi- 
plied by two to eliminate his Earth shielding factor. 

At first thought, an assumption might be made that the results of this analysis 
should have indicated a substantially greater influx than had previously been thought to 
exist. For, the analysis shows a much greater number of meteoroids of low velocity 
relative to those of higher velocity. The additional meteoroids at low velocity might be 
considered as additive to those that had previously been counted. However, by study of 
figure 5, it may be seen that such is not the case. By previous methods, only the two 
or three lowest plotted points would have appeared as lying within a straight line with 
a slope corresponding to a value of /3 = 1. 34. All the other plotted points would have 
fallen off below that line. By previous methods it was essentially assumed that such 
fall -off was due to failure of observation of faint meteors and the slope for the straight 
line was based on the lowest points. But as even the slowest meteoroids were observed 
for the largest masses, approximately the proper slope and the proper level of the line 
had been determined on the basis of those lowest points . 

Hence, the establishment of a substantially reduced average geocentric velocity of 
meteoroids in this analysis is actually in harmony with the previous concept as to mass 
distribution and influx rate, and it reinforces that concept by the demonstration that a 
straight-line relation actually does exist at least for the eight heaviest of the 10 mass 
groups analyzed. 


CONCLUDING REMARKS 

After completion of the analysis that has been presented herein, references 20 and 
21 came to the author’s attention. In reference 20, Erickson derived a velocity histo- 
gram that closely resembles the adjusted velocity histogram for all masses presented 
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herein, with an average velocity of 16. 5 kilometers per second (2.7 km/sec greater than 
the value derived herein) . His method was roughly complementary to that used herein, 
inasmuch as he divided the meteor sample into velocity groups rather than mass groups . 
In reference 21, Dohnanyi showed a similar "semi-empirical” distribution of velocity 
v TO , differing from that obtained herein for velocity v G only slightly more than 
Erickson’s distribution. Dohnanyi’s average velocity is 19.2 kilometers per second. 

The weighting factor that was developed in reference 2 is supported by the general 
agreement of a histogram for unaccelerated velocity of all meteors, adjusted by that 
weighting factor, with an unadjusted histogram for unaccelerated velocity of heavy 
meteors. That agreement, considered together with the partial agreement of the ad- 
justed histogram for all masses with the unadjusted histograms for 10 mass groups, 
supports the assumption of a single velocity distribution for all masses within the 
photographic range. Also, the analysis for cumulative mass distribution shows a linear 
relation on a full logarithmic plot. 

It is true that the assumption of a single -velocity distribution for all masses as well 
as the straight-line relation and the value of 0 = 1.34 were used in the derivation of the 
weighting factor in reference 2. In appendix E an argument is presented that each of 
these assumptions is confirmed independently of their use as assumptions in reference 2 . 
Hence, the results of the present analysis are of major value in establishing the depend- 
ability of the weighting factor. At the same time, the establishment of a basis for con- 
fidence in the weighting factor automatically establishes a basis for confidence in the 
average unaccelerated geocentric velocity of about 13. 8 kilometers per second obtained 
with its use. This average velocity (not greatly different from the averages of Dohnanyi 
and Erickson) is less than half the raw average of 30 kilometers per second, which has 
had much use. The greatly reduced average velocity should have an important effect on 
estimates of necessary armor thickness. Further study in this connection would seem to 
be well justified. 

As stated in reference 2, dissension exists as to whether luminous efficiency is pro- 
portional to meteoroid velocity relative to Earth’s atmosphere. The weighting factor 
^w(l) was derived on that basis. In appendix F an argument is presented that the present 
analysis tends to support the linear dependence of luminous efficiency on velocity. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 21, 1968, 

120-27-04-36-22. 


APPENDIX A 
SYMBOLS 


'cum(i) 


F(Zr ) 0 


f(v r ) 


effective capture cross section of Earth for meteoroids (larger than Earth’s 
cross section because of gravity focussing) 

constant used with varying significance as made clear in context 

constant used with varying significance as made clear in context 

count of meteors within mass group i, adjusted only for effect of gravitational 
focussing 

fully corrected cumulative count for influx of meteoroids of mass greater 

o 

than lower limit for mass group i as defined by eq. (20), particles m" 
sec"* 

fully corrected count of meteoroids within mass group i defined by eq. (18) 

tentatively corrected count of meteoroids within mass group as defined in 

2 1 

eq. (15), particles m" sec 

constant used with varying significance as made clear in the context; also, 
particle flux at infinite distance from Earth, as used in eqs. (Bl) and (B2) 

statistical function of zenith angle of meteor and position of meteor within 
field of view of cameras, explained in detail in ref. 2, nondimensional 

influx rate of meteoroids of mass greater than specified value m, particles 
m”2 sec"* 

frequency of occurrence of meteors of geocentric velocity Vq as defined in 
eq. (C3) 

subscript for the symbol equal to 1, 2, 3, or 4, and indicating that 

the symbol is defined by the equation of that number 

radius of a circle of area A, km 

radius of Earth, km 

velocity of escape from Earth at position where meteoroid impacts Earth 
atmosphere, km/sec 

velocity meteoroid would have had relative to Earth if Earth’s gravity had 
not affected that velocity, km/sec 



I ■“ 

v M velocity of meteoroid relative to Earth atmosphere at position of first infinitesi- 
mal encounter with atmosphere, km/sec 

hypothetical geocentric velocity of particle at infinite distance from Earth, 
km/sec 

x distance from y-y axis in pair of rectangular coordinates, used with varying 
significance as made clear in the context. 

y distance from x-x axis in pair of rectangular coordinates used with varying 
significance as made clear in the context 

Zj^ angle of path of meteor to zenith, deg or rad 

a factor of proportionality in eq. (24) 

or* factor of proportionality evaluated in eq. (11) 

/3 exponent of mass m in eq. (24) for influx rate of meteoroids 

<p w ( n ) weighting factor defined in equation having number indicated by subscript n, 
nondimensional 


*e> 

’’Jo 


revised correcting factor for counts of meteors within mass group i defined by 
eq. (17) 

correcting factor for counts of meteors within mass group i (identical with <p 
of appendix C) 
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APPENDIX B 



where Ig is number of particles intercepted by Earth, rg is radius of Earth and 
atmosphere, c is particle flux per unit area at infinity (in parallel paths), v is veloc- 

V 

ity of escape from Earth, and is particle velocity at infinity, primed here to avoid 
confusion with another significance assigned to v^ in the main text. With the aid of 
equation (Bl), an adjusting factor may readily be formulated for use in the main text to 
adjust an actual count of particles at any particular geocentric velocity, as observed 
within the atmosphere of Earth, to yield a value that would properly apply to a gravity - 
free Earth. Such an adjusted count would apply to a space radiator located approxi- 
mately in a position along the path of Earth in its orbit, but at a position far from Earth. 
Before such adjusting factor is formulated, however, a verification of the applicability 
of equation (Bl) will be undertaken. 

In the derivation of equation (Bl) in reference 7, the gravitational influence of the 
sun was ignored and particles were treated as emanating perpendicularly from a circu- 
lar area such as A within a plane surface located at an infinite distance from Earth 
(see sketch (a)). The reality, of course, is that particles approach the Earth in orbits 
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that are conic sections having foci at the center of the sun. With neglect of the sun’s 
influence, the ratio of the number of particles reaching Earth to the number that would 
have reached Earth without the effect of Earth gravity would be the ratio of the area A 
in sketch (a) to the projected area of the Earth and its atmosphere. But, if equation 
(Bl) is correct, then this same ratio should be equal to the left-hand side of equation 
(Bl). Hence, from equation (Bl), rewritten as 



(B2) 


a ratio of radii may be expressed as 



where r^ is the radius of the circular area A. 

Equation (B3) can easily be derived for the condition assumed, with reference to 
figure 6. As circular symmetry exists about an axis extending from Earth's center to 
the center of the circular area A, it is necessary to consider only a single particle, 
leaving some point within the perimeter of the area A in a direction parallel to the 
axis of symmetry. Such a particle, being marginal from the standpoint of collision with 
Earth's atmosphere, should just graze the atmosphere tangentially. It is well known 
that a particle at an infinite distance from and having a component of velocity toward an 
attractive center will approach that center along one branch of a hyperbola having that 
center as its focus. This condition is illustrated in the figure. The hyperbolic orbit 
is assumed to be symmetrical relative to an x axis with origin midway between the 
two branches of the hyperbola. 

The general equation for such a hyperbola is 

2 2 

I (B4) 

a 2 b 2 

where a and b are constants. It is well known that 

a 2 + b 2 = c 2 (B5) 


35 



(a) Sun’s gravitation neglected. 

Figure 6. - Effective capture area of Earth for slow particle. 


where c is the distance from origin to the focus of either branch of the hyperbola. 

The physical significance of a, b, and c, as shown by either of the two congruent right 
triangles in figure 6(a), is also well known (see any textbook in analytic geometry). From 
the figure it is seen that r^ is equal to b. The values of a and b may be determined 
from conditions that must exist at perigee, the point where the applicable branch of the 
hyperbola intersects the x axis. Those conditions are that (1) the marginal particle 
must reach a minimum distance from the center of Earth equal to r^, (2) the gravita- 
tional and centrifugal forces acting on the particle must extend in exactly opposite 
directions, and (3) because of item (2), the gravitational and centrifugal forces must be 
equal. 

A foundation has now been laid for derivation of equation (B3), the equivalent of 
equation (Bl). The derivation may be simply performed with use of well known principles 
of physics and analytic geometry. The details are omitted. 

Difficulty is encountered in a rigid demonstration that equation (Bl), (B2), or (B3) 
applies even approximately to the real condition involving the sun's influence. However, 
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it can be shown that application of these equations is reasonably consistent with a method 
of approximate solution of the 3-body problem that has given good results in practice. 

Knip and Zola (ref. 23) have shown that good approximate solutions are obtained by 
the following method for trajectories of small objects from positions near Earth to 
positions far from Earth or the reverse. Here the expression "near Earth" applies to 
positions such that the gravitational influence of the sun is negligible compared with that 
of Earth. The expression "far from Earth” applies to positions such that the gravita- 
tional influence of Earth is ngeligible compared with that of the sun. The method is 
composed of the following procedures: (1) a "sphere of influence" of Earth is defined 
as a moving sphere in space, always concentric with Earth's spherical surface, and 
having a radius equal to 145 Earth radii, (2) when the small object is within Earth's 
sphere of influence, it is treated as being affected by Earth’s gravitation alone, except 
that the motion of the object is defined relative to a coordinate system with origin at 
Earth center and moving with Earth center in orbit about the sun, (3) when the small 
object is outside Earth's sphere of influence, it is treated as being affected by the sun's 
gravitation only, and (4) when the small object passes into or out of Earth's sphere of 
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influence, the velocity and direction of motion are treated as unchanged relative to either 
Earth or sun. 

In application of the method described to the gravity focussing effect, it is necessary 
to consider particles entering Earth’s sphere of influence in parallel paths at a given 
velocity. The condition is illustrated in figure 6(b) . For the same value of v^, the 
hyperbolic orbit is exactly as in figure 6(a). However, a marginal particle, destined to 
graze the Earth atmosphere, will enter the Earth sphere of influence traveling in a 
direction tangent to the hyperbolic orbit at the distance of 145 Earth radii from Earth 
center. If thereafter it traveled an undeflected path, it would follow the tangent line and 
would cross over the x axis at a position to the right of the origin. The circular 
source area A must now be treated as at a distance of 145 Earth radii from Earth center 
rather than at an infinite distance. A particle emanating from its center, in a direction 
parallel to the tangent line, would follow a straight-line path toward Earth center. 

Hence, the radius of the source area A must now be the perpendicular distance of 
the tangent line from Earth center. With the newly defined value of r^, the concen- 
tration ratio due to the focussing effect of gravity would still be expressed by the left- 
hand side of equation (B3), but not by the right-hand side of that equation. 

No analytical expression has been found to substitute for the right-hand side of 
equation (B3). However, for any particular value of v^ a value of the concentration 
ratio (r^/r^)^ ma y rea ^ily calculated with use of elementary principles of analytic 
geometry. For example, with v'^ =1.0 kilometers per second the value of (r^/r^)^ 
may be calculated as 67.15. For the same value of v’^, equation (B2) gives a value 
124.21. Hence, it is seen that equation (Bl) or (B2) would provide too large a concen- 
tration ratio by approximately 85 percent for v^ = 1 kilometer per second, if the sphere- 
of- influence treatment were considered to be strictly correct. 

Computations and comparisons were made in the manner described for each integral 
value of v^ from 1 to 80 kilometers per second. The results, in percentage over- 
estimate by equation (Bl) or (B2), in comparison with the sphere-of-influence treatment, 
are plotted in figure 7. The percentage overestimate reduces to less than 4 percent at 
a velocity v'^ = 5 kilometers per second. But, in figure 1, it is seen that particles with 
values of v^ = 5 kilometers per second or less are only a small fraction of the total. 

From the foregoing treatment, it appears that use of equation (Bl) is well justified 
for the purposes of the analysis conducted in the main text. Accordingly, suitable 
adjusting factors <P W ( 2 ) an(1 ^*w(3) wil * now be formu l ated on the basis of that equation. 
The desired factors will operate on the flux within Earth atmosphere observed at a 
given value of v^, with and without adjustment for the photographic bias of meteoroid 
mass and velocity, to produce a flux as it would exist without gravitational focussing. 

If c represents the flux corrected for the effect of gravitational focussing and c Qb 
represents the observed flux, with or without the adjustment for mass and velocity, 
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then 


c ob = 


7rr A 2c 


ITT 


E 


(B6) 


From equations (B2) and (B6), 


"'ob 


— 1 + 



(B7) 


or 


c = 


"ob 



(B8) 


Now, in the main text, a weighting factor is expressed in equation (1), by 

which an actual count of particles at any particular velocity may be multiplied to yield 
a value corrected for mass and velocity bias, and a weighting factor ( P W ^ equal to 
unity is shown by equation (4), use of which in the computer program leaves the actual 
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count unaffected. If the actual count is designated c^, then application of the factor 

«Vd or <°w(4) y lelds 


c ob “ c ob^w(l) 


(B9) 


or 


c ob c ob^w(4) 

in which c ^ is a count that may be substituted into equation (B8) to yield a count c 
corrected for the gravity focussing effect. Revised correcting factors <P W ^) anc * 

<?w(3) will now be defined, such that 

■\ 

c = c ob^w(2) 

or > (BIO) 


c = c ob^w(3)J 


Then, from equations (B8) to (BIO), 


c 




^w(l) 




or 



(Bll) 


L J J 

Now, in the main text, Vq is equivalent to v^' as used in this appendix. If such 
substitution is made in equations (Bll), and if the values of and from 

equations (1) and (4) are also substituted, equations (Bll) express as c/c^ the 
weighting factor of equation (2) or </? w ^ of equation (3). 
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APPENDIX C 


DEVELOPMENT OF OFFSET LOG-NORMAL EQUATION 
FOR VELOCITY DISTRIBUTION OF METEOROIDS 


The curve in figure 2 is a parabola centered on the y axis as shown. The equation 
of that parabola is 


x = ±0.409V=T (Cl) 

where x and y are linear distances measured from the axes x-x and y-y as shown, 
in a unit equal to log 10 as plotted on the full logarithmic chart. It can readily be 
shown that this parabola represents a log-normal velocity distribution, offset 1. 5 kilo- 
meters per second toward higher velocities (ref 8) . 

Choice of a unit of length equal to log 10, for convenience, in the formulation of 
equation (C 1) and in the construction of the parabola in figure 2 commits equation (C 1) 
to sue only in conjunction with logarithms to the base 10. For it is only to the base 10 
that log 10 can have a value of unity. In the measurement of x and y in figure 2, the 
unit of length log^Q 10 is greater than the unit of length log g e by the factor log g 10, 
because in the first case the distance from 1 to 10 on the logarithmic scale is the unit of 
length, while in the second case the distance from 1 to e on the logarithmic scale is 
the unit of length. If x and y in equation (Cl) each have the dimension l , then the 
constant 0. 409 must have the dimension Hence, if the unit of l is changed by a 

factor (log g 10)”'*', the value of the constant 0. 409 must be multiplied by the one -half 
power of the reciprocal of that factor. Accordingly, equation (Cl) must be changed to 

x - ±0. 409 ^log e 10 (C2) 


Equation (C2), then, may be used in conjunction with natural logarithms. 

The standard equation for a normal distribution of velocity Vq would be 


f(v G ) = a* exp 


1 

2 



2 


(C3) 


where f (v G ) represents the frequency of occurrence of any velocity v G in terms of 
fraction of total population having velocity v Q ±t£ dv Q , divided by dv G , and where a, b, 
and c are constants. For a log-normal distribution in terms of v G (offset 1.5 km/sec) 
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equation (C3) becomes 


or 


or 


, % , lP°g e ^ V G + 1 * 5 ) “ lo S e b T 

f(vQ) = a. exp ^ 5- 

2 _ c 


(C4) 


, v i flog (v G + 1.5) - log b 

log e f( VG ) = log e a — 


(C5) 


lo S e ( v G + 1 - 5 ) " lo g e b = \^ c [ lo S e a " log e f ( v G^ ^ 


(C6) 


Now if the x-x axis is located at an ordinate value equal to log a and the y-y axis is 
located at an abscissa value of log b (both on the full logarithmic plot of fig. 2) and if 
the values of x and y are measured linearly from these axis, in units equal to log. e 

V 

on the scale of figure 2, then 


x = log e (v Q + 1.5) - log e b 


(C7) 


and 


y = log e f ( vq) - log e a (C8) 

With substitution of x and y from equations (C7) and (C8), equation (C6) becomes 

x = ±-yJT c-y^y (C9) 

By comparison of equation (C9) with equation (C2), it is seen that 

t/logp 10 

c = 0.409 I - = 0.44 (CIO) 

V5 

From the positions of the x and y axes in figure 2 and with use of a fractional value 
rather than percentage value for frequency, 

log e a = log e 0. 07 (Cll) 
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and 


log e b =log e 11.7 


(C12) 


Finally, equation (B13) may be written from equations (C4) and (CIO) to (C12). 
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APPENDIX D 


ADJUSTING FACTOR FOR COUNT OF METEORS WITHIN A MASS GROUP 


It will be demonstrated that a valid adjusting factor for a count of meteors within a 
particular mass range is obtained as the antilogarithm of the vertical shift necessary to 
bring the higher -velocity regions of an unadjusted velocity distribution for the mass 
group into agreement with the correct distribution, when both distributions are plotted 
on a semilogarithmic scale, as in the case of figure 3(e). 

Within a given mass group, let f Q (vj) be the observed fraction of the observed 
total count, and let f r (vj) be the correct fraction of the correct total count of meteors 
having velocities within the intervale v^ ± 6, where 5 is some small number. Assume 
v^ is high enough that no meteors within the mass group and within the velocity range 
v^ + 6 will go unobserved because of faint luminosity. Then 


where k Q is the observed number of meteors and k r is the correct number within the 
mass group and within the velocity interval v^ ± 6. 

Let A q be the total observed count and A r be the total correct count for the mass 
group. Then 


w-f 

A o 


(D2) 


f 


r 



P3) 


Now, the desired adjusting factor, by which A Q should be multiplied to obtain A f , is 



(D4) 


From equations (D2) to (D4), 


<P = 


£ q( v 1> k i 
f r < v l> k c 


(D5) 
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From equations (Dl) and (D5), 


<P = 


f o < v l> 
f r ( y i) 


(D6) 


or 


log <p = log f 0 (Vj) - log f r (v 1 ) (D7) 

For semilogarithmic plots, equation (D7) indicates that, for any velocity v^ where 
the observed count is correct, the needed adjusting factor will be the difference in level 
of the observed and correct curves at that velocity v^. Hence, if theoretical reason 
exists to believe that the observed count is correct throughout a particular velocity 
region, then the logarithm of the desired adjusting factor can be obtained by measuring 
the vertical shift required to bring the observed distribution curve into agreement with 
the correct curve throughout that velocity region. 

But good reason has been given to expect that the observed curves will be correct 
throughout the velocity region above some critical value at or above which failure to 
observe meteors would not occur for the mass group. Hence the method should be 
applicable in the comparisons used in the main text. 
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APPENDIX E 


EFFECTS OF ASSUMPTIONS REGARDING MASS AND VELOCITY DISTRIBUTIONS 

A question of circular logic exists in the indication that /3 = 1 . 34 in figure 5 and 
the indication of identical velocity distributions for different masses in figures 3(a) to (i) . 
Both these conditions were assumed in the development of the weighting factor 
of equation (1) (eq. (75) of ref. 2). 


Mass Distribution 

Consideration will now be given to the question whether an iterative procedure that 
could be performed would converge to a value of approximately 1.34 for /3. That pro- 
cedure would involve repetitions of the work shown in reference 2 (pp. 31 to 34) and the 
work presented in this paper between equations (1) and (24). In each repetition, the 
value of j3 obtained from the preceding repetition would be used as a starting value in 
reference 2 (p. 31). Comparisons of histograms would be made in a new figure 3, and 
a new figure 5 would be constructed. 

The suggested iterative procedure would be excessively long and laborious. 

However, the question of convergence to a value of j3 near 1. 34 can be examined without 
the need of actually performing the procedure, as follows. 

In description of some details of the iterative process, the factors (cos Z-r,)"®' 

0 730 

and F(Zp) a y in equation (1) will be omitted, because an adjusted histogram for all 
masses (not shown in this report) with the weighting factor according to the equation 

K ( i )= v ~' 22 w 

does not differ appreciably from the adjusted histogram for all masses obtained with 
<^>w(i) according to equation (1). 

In the analysis that has been reported here, in determining the value of j3, the 
author depended on physical measurement of slope from a carefully constructed full- 
logarithmic plot, which has been reproduced as figure 5. This method was preferred to 
a least-squares fit because it allowed judgment in weighting of the plotted point for mass 
group 8 with allowance for the possibility that the fall -off evident for mass groups 9 and 
10 may actually have started with group 8. But, for the purpose of examining the question 
of convergence, a rigid mathematical determination of slope will be needed even if it 
does not necessarily give the best result. For that purpose, the straight line providing 
a least- squares fit to the data will be used for the first eight mass groups. Groups 9 and 
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10 will be disregarded. For the data shown in table II and plotted in figure 5, such a 
least-squares fit was obtained with a slope of -1. 32, which is only 0. 02 different from 
that determined visually. 

Now, as an example, if a value of /3 = 1. 25 had been assumed between equations (59) 
and (60) in reference 2, then the exponent of v^ in equation (75) would have been -3.96. 
The lower absolute value of this exponent would have changed the shape of the adjusted 
velocity distribution for all masses used for the comparisons in figure 3. But the change 
in shape would not have been enough to prevent the comparisons, as will be explained. 
Also, the effect of the changed assumption regarding /3 on figure 5 will be estimated. 

In the comparisons, instead of the adjusted histogram for all masses, use could 
have been made of a smooth curve representing f(vQ) according to equation (5) (times 
100). Some analogous smooth curve would exist with use of a weighting factor like 
<p’ w (i) of equation (El) but with the exponent -3.96. Relative influx rates shown by 
that analogous smooth curve would be proportional to those given by equation (5), but 
multiplied by v 0^. 22-3. 96) or v 0.26 

Now, with mass group 5 as an example, it is seen that the agreement in shape of 
the two histograms in figure 3(e) occurs approximately within the region of values of 
Vq from 28 to 58 kilometers per second. (Hereinafter, for any mass group i, the range 
of Vq within which the shape of the histograms agree will be referred to as the matching 
range.) The analogous smooth curve for the assumption /3 = 1.25 would not fit the raw 
histogram in this figure within the matching range so well, but it could have been fitted 
approximately, with a cross over somewhere within the matching range (values of Vq 
from 28 to 58 km/sec). 

The percentage of total influx of meteoroids shown by the analogous smooth curve 
at any velocity Vq ±0. 5 kilometers per second must be given very closely by the 
equation 


f,(v G ) (%) - 



Equation (E2) indicates an influx of 0.888 and 0.0104 percent of the total with velocities 
Vq at 28 and 58 kilometers per second, respectively, which are higher on the loga- 
rithmic scale than the values 0.771 and 0.00757 percent given by equation (5) by dis- 
tances equal to log (1. 152) and log (1. 374). Hence, the analogous smooth curve for the 
assumed /3 = 1.25 would require a smaller shift to produce approximate agreement 
with the raw histogram within the matching range for mass group 5 by a distance approx - 
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imately equal to the average of log (1. 152) and log (1. 374), that is, by the distance log 
(1.258). Accordingly, the shift cp'^ in table I would be reduced from log (6.70) to log 
(6.70/1.258), or log (5.326). (Hereinafter for any mass group i the reduction of shift, 
like log (1.258) for i = 5, will be called the shift decrement.) 

Similarly, <p ^ and (p^ * n table I would become log (1. 32) and log (1. 14), the 
last term in equation (17) would become -log (1.223), the corrected shift in table n 

would become log (4. 36), and the fully corrected count Cg would become 1.72x10“*^. 

Reduced values of C i were also found in the same manner for mass groups 1 to 4 
and 6 to 8. The resulting fully corrected cumulative counts C cum ^ for groups 1 to 8, 
were 3.94, 7.22, 14.66, 27.46, 44.66, 86.81, 167.5, and 346. 5, all multiplied by 
Hr 15 . The straight line providing the least-squares fit to the logarithms of these data 
has a slope of -1. 27, only 0.05 different from the result of -1. 32 for the data in table n. 
But, as the new assumed value of £ s 1.25 is less than the old assumed value /3 = 1. 34 
by 0.09, it may be seen that convergence has occurred by 44. 5 percent of the difference 
between the two assumed values. 

The foregoing procedure was also performed with an assumption /3 = 1. 43. In this 
case, shift increments occurred instead of shift decrements. The resulting values of 
corrected cumulative counts for mass groups 1 to 8 were 3.94, 7.22, 16.07, 31.36, 

52. 33, 112. 3, 224.5, and 484.0, all multiplied by 10"*®. The straight line providing 
the least-squares fit in this case has a slope of -1. 38. Here, convergence occurred by 
33. 3 percent of the difference between the assumed /3 = 1.43 and the original value 
0 = 1.34. 

The three pairs of values (assumed /3 equal to 1.25, 1.34, and 1.43, with resulting 
indicated values by least squares of 1.27, 1.32, and 1.38) indicate that convergence 
would occur to a value of 1.30. With indicated values by the graphical method in each 
case higher by 0.02, the indicated values of /3 for the three cases would be 1.29, 1. 34, 
and 1.40, and the indicated value to which convergence would occur is very closely 1.34. 

It may now be seen that convergence occurs for the following reasons: 

(1) If the comparisons of histograms had involved the same matching range for all 
mass groups, then assumption of j3 = 1.25 (or any other value) would have decreased 
all values of log by the same shift decrement; consequently, both terms in the 
right-hand side of equation (17) would have been decreased by the same amount. Hence 
all values of <p ^ from equation (17) would have been unchanged and no change would 
have occurred in any value of C cum ^. Any change in the values of C cum ^ and 
consequent change in indicated value of 0 would have to be caused by progressive 
change in the matching range and consequent progressive change in the shift decrement 
from one mass group to another. 

(2) The increase in ordinate between mas groups i - 1 and i in figure 5 is equal 
to 
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(E3) 


A1 °S C cum(i-l,i) =1 °S 


1 + 


C i 


'cum (i-1) 


Any change in matching range between group (i - 1) and group i, with consequent change 
in shift decrement, may change the value of and hence the value A log C cum ^_j ^ . 
But note in table n that and C cum (i_i) are generally nearly equal. Hence, the 
values of C cuma . s) , C cum(1 . 3); and so on, are comparatively unimportant parts of 
Ccum(i-l) as ^ a PPears in equation (E3) . Consequently, if a change in shift decrement , 
occurs between groups (i-1) and i, then between groups i and (i+1) that changed shift 
decrement will apply almost as much to the denominator in equation (E3) as to the 
numerator. It follows that substantial changes in the value A log C cum ^_j ^ will 
occur only because of change in the shift decrement between groups (i-1) and i. The 
effect of that change in shift decrement will rapidly be lost for groups (i+1), (i+2), and 
so on. 

(3) As an example, the value of 0 indicated by the interval between mass groups 
4 and 5 only (from values in table n) is 


$ I°gl 0 (4.818xl0- 14 )-log 10 (2.916xl0- 14 ) _ 0 2181 _ i ^ ^ 

( 4 ’ 5 > log 1Q (0. 1025) - log 1Q (0. 1603) 0.1942 


Now, from equation (E4), and by the reasoning in items (1) and (2), any change in shift 


decrement AD 


(4,5) 


between mass groups 4 and 5 will change the value of 0 


(4,5) 


to 


0 


(4, 5) a 


0.2181 - AD 
0.1942 




(E5) 


or 


A/3 (4, 5) " 0(4, 5)a ' 0(4, 5) 


~ AP (4,5) 

0.1942 


(E6) 


The median values of v c within the matching range for groups 4 and 5 were 41 and 43 
kilometers per second. The corresponding values of v ro are 42. 5 and 44.4 kilometers 
per second. Then, as the denominator of equation (E2) is constant, the change in shift 
decrement, by equations (5) and (E2), is 
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iD (4, 5) - ^lO (^i)°' 2S - l^wd' 045 ) 0,26 


(E7) 


So, from equations (E6) and (E7) 


a ^ 5 ) = -0.0267 (E8) 

This result is about half the reduction 0.05 in the indicated value of /3 obtained by 
the least-squares fit for the first eight mass groups, probably because the change in 
median value of the matching range between mass groups 4 and 5 is somewhat less 
than in some of the other cases and because the conclusion in item (2) is only approx- 
imate . 

(4) By review of the foregoing and the pertinent parts of reference 2, it may be 
seen that the exponent 0.26 in numerator and denominator of equation (E2) is equal to 
the factor 2.87 (ref. 2, eq. (59)) multiplied by the reduction of the assumed value of 
/3 (from 1.34 to 1.25 in this case); that is, 

p = -2.87A/3 & (E9) 

where p is an exponent that might be substituted for 0.26 in equation (E2) and A/3 

d 

is a change in assumed value of j3 from the originally assumed value of 1.34. Hence, 
with any assumed value of /3, the change in shift decrement between the fourth and fifth 
groups as an example would be given by the following modification of equation (E7) 

AD (4) 5) = log 10 (1.045)- 2 - 87 A/3 a (E10) 

and from equations (E6) and (E10) the change of indicated /3 would be, approximately, 

, -2.87 A j8 

A P/a log 1Q (1.045) a = -0.28 A/3 (Ell) 

0.1942 iU a 

Equation (Ell) indicates that, whatever value of /3 had been assumed, for the interval 
between mass groups 4 and 5 convergence would have occurred about 72 percent of the 
way toward the value 1.32 as shown by the least-squares fit, or the value 1.34 obtained 
by the graphical method. 
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It appears now that the question of circularity is eliminated relative to the indicated 
value of j3. It must still be examined regarding the matter of the same velocity distri- 
bution for different masses. The indication that * 1. 34 is still dependent on 
correctness of the assumption regarding uniformity of velocity distribution. 


Velocity Distribution 

As discussed earlier, the comparison of histograms in figure 3 does not agree 
exactly with the assumption of a single velocity distribution, as it affects values of Vq 
below about 20 kilometers per second in comparison with higher values. However, the 
assumption is not contradicted, even as an exact condition. Discussion here will be 
understood to apply only to velocity ranges in which confirmation exists by the compar- 
isons in figure 3. 

As a starting point for discussion, let it be supposed that the condition assumed 
actually does exist, just as shown by figures 3(a) and (b), for the two heaviest mass 
groups, and by figures 3(c) to (i) for other mass groups throughout the velocity region 
above a critical velocity in each case. All the histograms would then be correct as they 
are, because both assumptions used in the derivation of the weighting factor are 

then satisfied. 

Next let it be assumed that the influx rate, for the heaviest mass group only, is 
increased by an arbitrary percentage, but only throughout an arbitrary small velocity 
range from v^^ to Vq^)’ The level of the raw histogram of mass group 1 will be 
raised, within the region from to v^g) > by the same percentage as the increase 

influx. But, regardless of the value of 0, unless the range v q(i) to v q( 2) Resina 
region where the only incidence of meteors was within mass group 1, the level of the 
adjusted histogram for all masses, within that velocity region, will be raised by a 
smaller percentage. Hence, clearly the two histograms could not agree in shape after 
the change was made, even though the originally assumed condition would then be 
incorrect. 

The same reasoning would apply to any part of the histogram for any mass group. 

It is apparent, therefore, with only minor exceptions, that change in the velocity dis- 
tribution for any mass group would destroy the agreement of the histograms. Once the 
agreement was destroyed, it could be restored by changing the distributions for other 
mass groups, but only by changing in the same direction, that is, in the direction toward 
greater uniformity of the velocity distributions. 

This reasoning is not dependent on the manner of derivation of , except in the 

primary assertion that the comparison of histograms would be as it is if the velocity 
distributions were the same for all mass groups. 
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APPENDIX F 


RELATION OF LUMINOUS EFFICIENCY TO METEOROID VELOCITY 

In reference 2 a relation was used 


r = r 0 v (FI) 

in which r is luminous efficiency, v is instantaneous meteoroid velocity, and Tq is a 
constant. This assumption was involved in the derivation of equation (40) of reference 2, 
expressing a criterion for marginal exposure from the standpoint of discovery of a meteor 
on a photographic plate. Hence, it was involved in the derivation of <p of reference 2 
(eq. (75)). As stated in reference 2, equation (FI), with Tq a constant, is by no means 
universally accepted. For that reason, the mutual relation between equation (FI), with 
constant r Q , and the present analysis will be discussed. 

At least approximately Tq appears to be equal to a constant times v n , where n is a 
positive or negative constant or is equal to zero as treated in reference 2. If this condi- 
tion were substantially incorrect, then figure 4 of reference 2, for v M as the sorted pa- 
rameter, should have shown systematic variation of ordinate for various values of abscissa. 

Now, let it be assumed that n is not actually equal to zero. Then, if Tq is 
redefined as a different constant than before, equation (FI) becomes 

t = r^ 1 (F2) 


Then equation (40) of reference 2 becomes 

C marg = (m;/- 020 v 8l,42+n (cos Z R ) 1 - 080 F(Z B r I (F3) 

Now, however, m^ of equation (F3) cannot be the mass determined in reference 6. 
Relatively, without regard to absolute values, it must instead be 


m 


T 

OO 


m 


oo 


V 


-n 


(F4) 


where m^ is the mass shown in reference 6. Hence, in terms of m^, equation (F 3) 
becomes 
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C„„„_ = (rn 
marg x ° 


-n) 


1.020 


. 2 . 842+n 


(cos Z R ) 


1.060 


F(Z r ) 


-1 


= m l-020 y 2. 842-0. 02 n (cog Zr) 1.060 F (Z r ) _1 (F5) 

which, if the absolute value of n is not greater than 1 or 2, is not significantly different 
from equation (40) of reference 2 . 

Now C ... is a theoretically derived criterion for discovery of a meteor trace on 
marg 

a super Schmidt photographic plate. The ratio of the exponent of v^ to the exponent of 
m oo was confirmed closely in equation (45) of reference 2 by an empirical method. That 
method was based on 100 test meteors chosen in a manner not affected by in or m’ . 
The method is as good for one value of n as another in equation (F4) . But not so with 
the weighting factor <p w of equation (75), which was developed in reference 2 (pp. 27 to 
34) with use of equation (45) and other considerations. Of the other considerations, those 
pertinent here are 

(1) the cumulative mass distribution 




(F6) 


(2) the assumption that the velocity distribution ffv^ (or f(vg)) is the same for all 
values of m^. 

Now it will appear later that equation (F6) may be used, with the same value of /?, 
with m^ or m^. But validity of the assumption of the same velocity distribution for all 
values of m^, under the assumption that m’^ is correct, needs some detailed consid- 
eration. For that purpose, alternative assumptions will be made as follows 

(1) fCv^) is independent of m^. 

(2) the cumulative frequency for m^, Fv >m j o , for a particular value of is 
independent of that value of v w . 

(Assumption (2) is a corollary of assumption (1)). Now the assumption actually used in 
reference 2, that ffv^) is independent of m^, does not necessarily follow from the 
assumption that ffv^) is independent of m^. The fact that it does follow will have to be 
demonstrated. 

Now assume any given value of m^ such as m^ and any two values of v^ such as 
Vj and Vg. Expressions may be written for two influx rates, 

N(m 1? v t ) = f(v t ) F >m ^ = 0!m 1 ~ /3 f(v 1 ) (F7) 
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and 


N(m 1 , v 2 ) = f(v 2 ) F >m ^ = am-f f(v 2 ) (F8; 

where N(m^,v^) is the influx rate for particles of mass greater than and 
having velocity equal to v±0.5 and N(m^,v 2 ) is the same function for v 2 . Let 

Nfm^vJ f(v 2 ) 

r’ = — tL = — (F9j 

Nfm^Vj) f(v x ) 

Now consider a similar comparison in which m^ is kept constant instead of m^. 
In this comparison, let the first influx rate be N(mj, v^) according to equation (F7) just 
as before. Then, for N(mj,v^), the value of according to equation (F4) will be 


"Ml) = m l V l 


(F10) 


For the second influx rate in this comparison let m^ have a value m 2 such that the 
corresponding value of m^ will be identical with m^^ of equation (F10), and let the 
velocity have the same value v 2 as in the eariler comparison. Then, from equations (F4) 
and (F10), 


m 2 = m °°(l) v 2 =m l 


(Fll) 


N (m 2 , v 2 ) = f ( v 2 ) F >m 9 = “ m 2 /3f(v 2 ) = 

* V2 


(FI 2) 


Hence, the new ratio of influx rates is 


N(m 2 ,v 2 ) /v 1 V p f(v 2 ) 
NOn^Vj) lv 2 y f(vj) 


(F13) 
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Finally, 



Equation (F14) is free of both and m^. It indicates that the ratio of influx 
rates at any two velocities is different with m^ constant than with constant, but 
nevertheless unchanging from one constant value of m^ to another. 

Hence, the assumptions used in reference 2 (pp. 27 to 34) are valid, and equation (75) 
gives a theoretically correct expression for <p w , even if the values of m^ are correct. 

It is only necessary to remember that the expression was derived on assumptions of uni- 
form velocity distribution for all values of m^ and a value of /3 = 1.34 applying to the 
cumulative distribution of m^. Hence, the validity of the weighting factor depends on 
confirmation of both those assumptions. 

But it is argued in appendix E that both those assumptions are confirmed by the 
present analysis. Hence, within the accuracy of the method, it seems that is a 

correct weighting factor, even under the assumption that n in equation (F4) is not 
equal to zero and under the corollary assumption that the correct mass is m^ of 
equation (F4). 

As far as has yet been shown, however, the value of /3 from figure 5 applies only 
to the values of m^. That it must apply also to the values of m^ may be seen as 
follows . 

If the values m^ are assumed correct, then the abscissa scale in figure 5 becomes 
logCm’^v 11 ). Now consider arranging a parallel abscissa scale as log m^. Consider in 
particular a value m^ = m^. As values of did not go entirely to zero, any value of 

= m^ in figure 5 must have lain within a region on the abscissa scale entending from 

log m^ = log(m lV ^ in ) (F15) 

to 

log m^ = log(m 1 vJJ 1 , ax ) (F16) 

where v min and v max are the minimum and maximum values of v^ actually encount- 
ered. The range of this region on the abscissa scale then is 
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R( m l> - 1 °g( m l'Cax - log ( m l'mln> = 1 o S P 222 ] ( p l 7 ) 

\ miiy 

It is seen that R(mj) is independent of the value of m^. 

Now consider the distance between similar ends of the range for two values of 
such as and mg. This distance is 


A log m^ = log(m 2 v^. n ) - log(m lV ^ in ) = log) 



(F18) 


But for the same two values of m^ and again with the minimum value of v , the values 
of m^ according to equation (F4) are 


m »(l) - “An 


(F19) 


m / 9 v = m«v m . 
°°(^) 2 mm 


n 


(F20) 


the same as before, and the distance between these two values on the abscissa scale is 


n n / U1 2 ’ 

A log m^ = log(m 2 v min ) - log(m 1 v m . n ) = logf — 


(F21) 


From equations (F18) and (F21), it is seen that the distance on the abscissa scale 
between m^ and m 2 is the same as the distance between m oo(i) anc * m oo(2)‘ 
scales are linear. 

So the abscissa scale for m^ (fig. 5) is the same as for m^, except that an 
uncertainty always exists as to just where to put a value of log within a fixed range 
of values of log m^. The same uncertainty exists as to the ordinate, within the range 
of values shown in the figure for logfm^v j n ) and log(m^v max ) . 

But, now, if the plot in figure 5 is assumed to extend from = 0 to m^ = 
the limited uncertainty as to where to place each value of abscissa and ordinate becomes 
of no consequence relative to the slope of the line that will be produced. Note that this 
theoretical device is not vitiated by the fact that the value of /3 will not remain fixed 
outside the mass range actually represented in figure 5. Hence the value of /3 must 
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equal 1. 34 for either or 

Now, returning to reference 2, let it be supposed that equation (F2) had been used, 
with n not equal to zero . Suppose further that equation (F3) had been used instead of 
equation (40) . Suppose that the mass values published in reference 6 had been the values 
of instead of m^. 

Then the empirical work between equations (40) and (45) should have confirmed 
equation (F3) , just as it actually did confirm equation (40) . But the work between 
equations (45) and (59) would have produced 


<P 


cor 


= F ( V 'oo Z R 


-/3[(2. 842+n)/l. 020] 


(F22) 


Then, the work between equations (59) and (75) could have produced even approximately 
the same weighting factor as before only if a value of /3 very different from 1. 34 
had been used. As both the weighting factor and the value of /3 have already been 
confirmed (even for n not equal zero), it seems to follow that a correct weighting 
factor could not have been found with use of equation (F22). 

At first thought, it may seem that both methods should be equally valid and that 
each should have its chance. Just as the value that was actually used is as valid 

for one condition as the other, then also equation (F22) and the replacement of equa- 
tion (75) in reference 2 that would have been derived from it should also be valid for 
both cases. 

But the new equation (75) would be valid for either case only if an analysis like that 
performed here confirmed, simultanteously, both the actual value of /3 used and the 
independence of velocity distribution on mass. As the present analysis did very 
closely confirm both those conditions, and did so simultanteously, it seems necessary 
to predict that use of a new equation for based on equation (F22), could not 

simultaneously confirm both conditions and that such a new weighting factor would 
therefore not be valid. 

Hence, within the accuracy of the method, the present analysis seems to support 
equation (FI) with constant r Q . Slight departures from this condition, particularly if 
they involved only a gradual change in t q approximately within the range of v^ below 
about 20 kilometers per second, would be possible. But, if later study shows that some 
low-velocity meteors were undiscovered within the, heaviest mass groups, then the 
possibility of error in the assumption that t q is constant will become even more 
restricted. 
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